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Abstract 

In the previous two parts of this series of papers, we introduced and studied a large 
class of analytic difference operators admitting reflectionless eigenf unctions, focusing 
on algebraic and function-theoretic features in the first part, and on connections with 
solitons in the second one. In this third part we study our difference operators from 
a quantum mechanical viewpoint. We show in particular that for an arbitrary dif- 
ference operator A from a certain subclass, the reflectionless A-eigenfunctions can be 
used to construct an unbounded self-adjoint reflectionless operator A on L 2 (M.,dx), 
whose action on a suitable core coincides with that of A. 

1 Introduction 

In this paper we study various quantum mechanical features of a large class of analytic 
difference operators that admit reflectionless eigenfunctions. Our analytic difference ope- 
rators (from now on AAOs) are given by 

A = exp(-id x ) + V a (x) exp(id x ) + V b (x), (1.1) 

where V a (x) and V b (x) are meromorphic functions with asymptotics 

lim V a (x) = 1, lim V b {x) = 0. (1.2) 

I Re ir I — *oo I Rex I ^00 

The notion 'reflectionless eigenfunction' refers to meromorphic functions W(x,p) satisfying 
the eigenvalue equation 

(AW)(x,p) = (e p + e~ p ) W(x,p), p £ C, (1.3) 

and having asymptotics 

I e lxp , Rex — > 00, 

W(x,P)~{ f lxp R 1.4 
I a(p)e p , Rex — > -co. 

We make extensive use of previous results in this series of papers, denoting Refs. [|J 
and pj by Part I and Part II, resp. In Part I we presented and studied a huge class 
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of reflectionless AAOs, but here we are concerned with a far smaller class. Indeed, our 
main interest in this paper is in associating with the AAO A a well-defined self-adjoint 
operator A on the Hilbert space 7i x , where we use the notation 

H y = L 2 (R,dy), (1.5) 

and we are only able to do so by imposing drastic restrictions on the spectral data in terms 
of which the coefficient functions ('potentials') V a and V b are defined. 

It should be mentioned at the outset that we are dealing with exotic territory. In 
Ref. H we studied in great detail some quite special reflectionless AAOs arising in the 
context of (reduced, 2-particle) relativistic Calogero-Moser systems. (See our contribu- 
tion Ref. Q to the NEEDS 2000 Proceedings for the unitary similarity transformation 
connecting Ref. || to the present framework.) From the findings reported in Ref. |3| it 
is already clear that a complete Hilbert space theory for reflectionless AAOs is not go- 
ing to amount to a straightforward extension of the well-known results for reflectionless 
self-adjoint Schrodinger and Jacobi operators. In Sections 3 and 4 of Part II we have 
summarized the latter results, and we have delineated restrictions on the spectral data for 
our AAOs such that their reflectionless eigenfunctions can be tied in with the Schrodinger 
and Jacobi counterparts. 

Roughly speaking, we impose similar restrictions in the present paper. We shall be 
quite precise in Section 2, but in this introduction we try and outline our results with 
a minimum of technical detail, as this might obscure the basically simple plan of this 
paper. Before sketching the latter, we add some general remarks yielding more context. 
To begin with, since we aim to associate with A a self-adjoint operator A on Ti x , it is 
natural to restrict V a and V b such that A is at least formally self-adjoint. Thus V b (x) 
should be real-valued for real x, and V a (x) exp(id x ) should be equal to its formal adjoint, 

[V a (x) exp(id x )]* = exp(id x )V a (x) = V a (x - i) exp(id x ), iGl. (1.6) 

Hence we need 

V b *(x) = V b (x), V a *(x)=V a (x-i), xeC, (1.7) 

where the * denotes the conjugate meromorphic function, 

f*(x) = m, xeC. (1.8) 

From now on we restrict attention to potentials satisfying (|l.7| ). Then a natural strategy 
would be to try and find a dense subspace C in TL X on which A is well defined and 
symmetric. Thus, C should consist of square-integrable functions f(x), x € M, that are 
restrictions to M of functions that have suitable analyticity properties for |Ims| < 1, so 
that there is an unambiguous meaning for f(x±i); then the function 

(Af)(x) = f(x-i) + V a (x)f(x + i) + V b (x)f(x), x£K, (1.9) 

should be square-integrable, and one should have 



(f,Ag) = (Af,g), f,geC. 



(1.10) 
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Assuming such a dense subspace has been isolated, one can try and study the existence 
and uniqueness of self-adjoint extensions. Indeed, the symmetric operator A on C is 
unbounded (due to the shifts), so it might not have any self-adjoint extensions or a (finite- 
or infinite-dimensional) family of self-adjoint extensions. 

In any event, assuming some self-adjoint operator A has been associated with A via 
this procedure, one can define its being 'reflectionless' solely in terms of time-dependent 
Hilbert space scattering theory, as follows. 

First of all, there is a natural 'free' dynamics exp(— UAq) with which the 'interacting' 
dynamics exp(— itA) can be compared. Indeed, the AAO 

A = exp(-id x ) + exp(id x ) (1-11) 

gives rise to an obvious self-adjoint operator Aq on Ti x , namely, the transform 

i = F Q MJ^ X (1.12) 

of the self-adjoint multiplication operator 

(M/)(p) = 2cosh(p)/(p), feV(M), (1.13) 

with maximal domain D(M) under Fourier transformation 

/oo 
dpe ixp f(p). (1.14) 
-oo 



(Recall our notation Ql.5| ).) 

Now assume that the (strong) limits of the operator family exp(itA) exp(— UAq) for 
t — > ±oo exist and have equal range. Denoting these isometric wave operators by W±, the 
corresponding 5-operator 

S x = W^W- (1.15) 
is unitary. Since it commutes with the free evolution exp(— UAq), its transform 

S P = T^S^q (1.16) 

to Tip is of the form 

(S p f)(p) = T(p)f(p) + R(p)f(- P ), f G Hp, (1.17) 

for certain functions T(p), R(p). Then the dynamics A is by definition reflectionless 
when R(p) vanishes identically. 

Our summary of these notions from time-dependent scattering theory (about which 
a wealth of pertinent information can be found in Ref. ||) serves a twofold purpose. 
First, it has enabled us to sketch a general scenario in which the concept of 'reflectionless 
self-adjoint AAO' makes sense and can be studied. Second, we actually follow a quite 
different strategy in this paper, but time-dependent scattering theory does play a crucial 
role. Thus we are now better prepared to sketch our special setting, and compare it to the 
above approach. 

The main difference consists in our definition of the self-adjoint operator A: It hinges 
on using the quite special ^4-eigenfunctions W(x,p). (Note that in the general setting just 
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sketched, eigenfunctions of the AAO A need not be and are not mentioned.) Specifically, 
the eigenfunction transform 

/oo 
dpW(x,p)f(p) (1.18) 
-oo 

plays a decisive role in defining A. 

We have already seen the simplest example of this approach. Indeed, for the AAO Aq 



( 1.11 ) we defined the associated Hilbert space operator Aq by using the ^o-eigenfunctions 
ex.p(ixp), cf. ( |1.14j ). The unitarity of Tq is crucial here: Invertibility of Tq would 

not be enough for ( |1.12| ) to give rise to a self-adjoint operator Aq. To compare with the 
general strategy, we mention that a domain Cq of essential self-adjointness as considered 
above is for instance given by J-q(C^°(M)). The point is, however, that the latter domain 
cannot be readily described in terms of the position space Ti x . Moreover, even for Aq 
there exists an infinite-dimensional family of distinct domains of essential self-adjointness 
yielding distinct reflectionless self-adjoint operators on Ti x . (This can already be concluded 
from the special cases studied in Ref. M, cf. also Ref. Q. The present more general case 
yields a much larger family, as shown at the end of Section 4.) 

To appreciate the latter state of affairs, and, accordingly, the choice involved in tak- 
ing (1.18) as a starting point, a crucial feature of ^4-eigenfunctions should be recalled: 



They remain eigenfunctions with the same eigenvalue after they are multiplied by an arbi- 
trary meromorphic function with period i. In particular, this entails that when an AAO 



of the form (Ll)-(1 1 2) admits a reflectionless eigenfunction W(x,p) satisfying ( |1.3|) -(|T 



it also admits a reflectionless eigenfunction W(x,p) with any other function d(p) in its 



asymptotics (|L4|). Indeed, we need only set 

W(x,p) = fi(x,p)W(x,p), (1.19) 

with 

fi(x,p) = (e 2 ™ + ~a{p)a{p)- 1 e- 2 ^ x ) / (e 2 ™ + e" 2 ™) , (1.20) 

to obtain a new eigenfunction with these features. 

Now there is no reason to expect that when the operator (|1.18j ) is unitary (or at least 
isometric) for a particular choice of W(x,p), it is still unitary/isometric for eigenfunctions 
W(x,p) as just described. Indeed, in the case of J-q it can be proved that multipliers of 
the form ( |l,20j ) destroy unitarity. But as already alluded to, in this case there does exist 
an infinite-dimensional family of i-periodic multipliers for which unitarity is preserved. 

In our approach, then, the Hilbert space features of the eigenfunction transform JF (1.18) 



are of primary importance. We are able to establish the relevant functional-analytic fea- 
tures by using the considerable amount of explicit algebraic and function-theoretic infor- 
mation gathered in Parts I and II. In particular, the surprising connection to classical 
iV-particle relativistic Calogero-Moser systems established in II Section 5 is instrumental 
in obtaining important additional information of the same character, whose derivation we 
have relegated to Appendix A. 

In outline, we solve the pertinent Hilbert space problems as follows. First of all, 
we choose the spectral data in terms of which W(x,p) is defined such that the trans- 



form T (1.18[) is a bounded operator on Ti p . This is already the case whenever W(x,p) 



has no poles for real x, which is a weak restriction. This choice also ensures that no 
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nontrivial Cg°(]R)-function is annihilated by T, cf. Lemma 21, As a consequence, we are 
entitled to define an operator A on the subspace 

V = TC^(R), (1.21) 

by setting 

ATf = TMf, /eC °°(K), (1.22) 
where M is defined by (1.13). (The relation to the AAO A is also clarified in Lemma 2.1 



A far more drastic restriction on the spectral data now ensures that W(x, p) has no poles 
for Imx £ [—1,0]. Note that this striking feature is generically destroyed when W(x,p) is 
multiplied by i-periodic multipliers fj,(x,p) with constant limits for |Rex| — > oo. (Indeed, 
by Liouville's theorem the latter must have poles in a period strip to be nonconstant. 
On the other hand, these poles might occur at the same locations as zeros (counting 
multiplicities) of W(x,p) in the strip Imx € [—1,0], in which case /j,(x , p)W (x , p) would 
still be pole-free in this strip.) 

Due to the absence of these critical poles, we are able to show that the operator A is 
symmetric on V . This involves considerable work, whereas the next step is quite easy: 
An application of Nelson's analytic vector theorem yields essential self-adjointness of A 



on V, cf. Lemma 2.2. Denoting the self-adjoint extension by the same symbol, we obtain 
a unitary one-parameter group exp(— itA) on the closure V of the subspace V. In general, 
this is a proper subspace of 7t x (that is, in general T is not onto Tlx), an d we now extend A 
provisionally to a self-adjoint operator acting in TL X by putting it equal to an arbitrarily 
chosen self-adjoint operator on the orthogonal complement . (At this stage we do 
not yet know that the latter space is spanned by eigenfunctions of the AAO A, so this 
provisional extension cannot be avoided.) 

Our next goal consists in handling the time-dependent scattering theory of the inter- 
acting dynamics exp(— itA), as compared to the free dynamics exp(— UAq). We do this 
in Section 3, the most important result being that the wave operators can be written in 



terms of T, cf. Theorem 3.2. From our explicit formulas it is then clear by inspection 



that T is an isometry. Moreover, they show that the 5-matrix S p ( 1.16| ) is the one ex- 



pected from time- independent scattering theory. (That is, the S'-matrix expected from 



the asymptotics (1.4) of the eigenfunction.) 

In Section 4 we complete our analysis by clarifying the state of affairs on This 
space is spanned by finitely many pairwise orthogonal eigenfunctions W(x,r n ), r n € tt), 
n = 1, . . . , N + , with distinct real eigenvalues 2cosh(r n ). Thus the definition of A can be 
completed by requiring that its action on V 1 - be equal to that of A, just as its action on V. 
The key to understanding is an explicit formula for TT* , which we obtain along the 
same lines as similar formulas for the special cases we studied in Ref. Q. (Since we have 
no duality properties available in the present general framework, we cannot proceed in 
this way to obtain the isometry formula T*T = 1. Instead, we exploit the isometry of 
wave operators, cf. Section 3.) 

We conclude Section 4 with an appraisal of some special cases. Of particular interest 
is the subclass for which no point spectrum occurs in the spectral resolution of A (cor- 
responding to N + = 0). This infinite-dimensional family has no analog for reflectionless 
self-adjoint Schrodinger and Jacobi operators. We use it to illustrate the ambiguity issue 
discussed above. 
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2 Essential self-adjointness on the domain "P 

We begin by recalling how our class of AAOs A Ql.l| ) and the associated reflectionless eigen- 
functions W(x,p) are obtained from 'spectral data' (r, /J,(x)). The vector r = (ri, . . . , rpj), 
with N € N*, consists of complex numbers satisfying 

e rm^ e ±rn f l<m<n<N, (2.1) 

and 

t c / (0,O, n = l,...,iV + , 

Imr n G < (2.2) 

\ (-vr,0), n = iV-iV_ + l,...,jV, V ; 

where 

N+,N- G {0,...,JV}, iV+ + iV_ = iV. (2.3) 
The vector = (pi(x), . . . ,/j,n(x)) consists of meromorphic functions satisfying 

fj, n (x + i) = fj, n (x), lim Hn{x) = c n , c n EC*, n=l,...,N. (2.4) 

\Rex\^oo 

These 'minimal' restrictions on (r, fi) are in force throughout this paper. (When the need 
arises, we specify additional restrictions.) 
Now we define a Cauchy matrix 

C mn = — — - — — , m,n = l,...,N, (2.5) 

g' m g ' n 

and a diagonal matrix 

D(x) = diag(d(ri,/ii;x),... ,d(r N , /ijvja;)), (2.6) 

where 

f#"^ ImpG(0,vr), 

y(a;)e -2i(p + »)x j l mp€( _^ ). 1 j 

Then the potentials 14, VJ, and wave function W are defined via the solution to the system 

(D(x) + C)R(x) = (, ( = (1, (2.8) 



by 



V a (x) = (l + e rn Rn(x)^J (l + eTnR ^ x + 0^ > ( 2 -9) 
Vb(x) = ^(i? n (x-*) -i? n (x)), (2.10) 



71=1 

/ N 



M^-e-^l-g^P^J. (2.H) 
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Of course, it is far from obvious that these definitions entail the eigenvalue equa- 
tion (03), but this is shown in I Theorem 2.3. In contrast, the asymptotics 



lim R(x) = 0, 

Re x— >oo 



lim R(x) 

Re x — ► — oo 



(2.12) 



easily follows from ([2.5|)- (|2.8| ), cf. also I Lemma 2.1. Using ( 2,12| ), one obtains the asymp- 



totics (1.2) and (O), with 



A' 



a{p) n 



n=l 



eP 



(2.13) 



cf. I Theorem 2.3. 

We are not able to associate a self-adjoint operator on 7i x to the AAO A unless we 
impose further restrictions on the data (r,[i). But to prove boundedness of the eigen- 
function transform T ( |1.18| ) and a few more salient features, we only need a quite weak 
assumption, as detailed in the next lemma. (Indeed, for generic spectral data satisfying 
( |2.l| )-( ^4|) the meromorphic function R(x) has no poles on the real axis.) 



Lemma 2.1. Assume that the solution R(x 
the operator T ( \1. 18\ ) is bounded. For all 
a meromorphic function satisfying 



to (2.8) has no poles for real x. Then 
the function (JF0)(:e) extends to 



E C °° 



A(F<j>)(x) = (TMfa(x) 



x E 



(2.14) 



where A is the AAO (LJ_) and M the multiplication operator ( LIS ). Moreover, we have 

Ker (T) n C^(R) = {0}. (2.15) 



Proof. In view of the asymptotics ( 2.12j ) and absence of real poles, the function R(x) is 
bounded for real x. Due to the restrictio ns (|2.2| ), the functions (e p — e~ Tn )~ 1 , n = 1, . . . , N, 
are bounded for real p. Hence W(x,p) (|2.1l|) is bounded for real x, p. Choosing 4>(p) E 
, we may write 



°0 



N 



n=l 



with 



(2.16) 



(2.17) 



cf. ( Oj ). Since Fourier transformation Tq is a bounded operator, and the multiplication 
operators occurring here are bounded, too, boundedness of J- follows. 

To prove the second assertion, we recall the easily verified fact that the Fourier trans- 
-function extends to an entire function. Since we have <f>, fa, . . . , <f>^ E 
and since Ri(x), . . . , Rn(x) are meromorphic, it is clear from (2.16) that (Tfa^) 
extends to a meromorphic function. The action of A on this function yields the meromor- 



form of a 

c n 



188 



S N M Ruijsenaars 



phic function 



N 



(Fofyix -i) -} y R n (x - i)(T o 0n)( 



x — i 



n=l 



+ V a (x) 



N 



(Fq4>){x + i) - y~] R n (x + i)(TQ(j) n )(x + 



+ VJ,(x) 



n=l 



(W)(a:)-13i2n(a;)(^b0)(a 



71=1 



(2.18) 



For all x for which the functions R(x), R(x ± £), 14 (a?) and T4(x) have no poles, this can 
be rewritten as the absolutely convergent integral 



(2vr)- 1/2 / d P [W{x -i,p) + V a {x)W{x + i,p) + V b (x)W(x,p)]<f>(p). 



(2.19) 



Thanks to the eigenvalue equation the function in square brackets amounts to 

2 cosh(p)W(x,p), yielding (2.14). 

To prove ( [Up , we assume <f> G Cg°(R) satisfies ^ = 0. By (|1|) we then have 



N 



{Fq4>){x) = y^ j R n {x){J : Q(t) n )(x). 



(2.20) 



n=l 



Consider now the function e ax R n (x), x G R, a > 0. Due to flgTp , it is bounded for 
x — > — oo. The Rex — > oo asymptotics of i?(x) can be sharpened to an exponential decay, 
so that e ax R(x) is also bounded at oo, provided a G [0,c], with c small enough. (The 
pertinent asymptotic decay easily follows from (|2.8j ), cf. I(2.41)-(2.42).) Now (.Fo</> n )(x) 
is a Schwartz space function, so it readily follows that the functions 



f n (z) = (27T)- 1 / 2 / dxe-^RnWiTo^ix), n = 1 



5"** J ! 



(2.21) 



are well defined for Imz G [0,c] and analytic for Imz G (0,c). Moreover, a dominated 
convergence argument yields 



lim f n (p + ia) = f n (p), 

uniformly for real p. 

From fl2~2(i| ) we deduce that the function 



n = l,...,N, 



(2.22) 



f(z) = (2vr)^/ 2 / dxe-^^ix) = £ f n (z) 

J -°° n=l 

is also well defined for Imz G [0, c] and analytic for Imz G (0,c). Since it satisfies 

lim f(p + ia) = 4>{p), 



(2.23) 



(2.24) 



uniformly for p G R, and 0(p) G Cq°(R) vanishes on an open set, Painleve's lemma entails 
(f> = 0. Hence (§T5|) follows. ■ 
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In view of ( |2.15| ), any vector in the subspace V ( |l.2l|) can be written as T§, <f> G C ( °° 



in a unique way. Therefore, the operator A ( 1.22j ) is a well-defined linear operator on V, 



whose action coincides with that of the AAO A, cf. ( |2.14 ). Obviously, A leaves V invariant. 
Moreover, we clearly have 

||i n ^|| < c"||^|| ll^ll, ra £ N, (2.25) 

with c > depending on supp(0). Thus V consists of analytic vectors for A. In view of 
Nelson's analytic vector theorem Q it now suffices for essential self-adjointness of A on V 
that A is symmetric on V . 

Next, we detail assumptions on the spectral data that suffice to prove this critical 
symmetry property. We do this in four steps, each of which adds a restriction. This 
enables us to use the less restrictive intermediate assumptions whenever we can show 
their sufficiency for the result at hand. (In most cases, however, we do not know to what 
extent these assumptions are necessary.) 

Our first step consists in imposing formal self-adjointness (|l.7j ) . This property is ensured 
by requiring that r\,...,r^ be purely imaginary and that the functions ie~ Tn fi n (x) be 
real- valued for n = 1,... ,N and real x, cf. I Theorem D.l. Our second step consists in 
requiring that /i(x) be constant. Thus our second assumption comes down to 

ir n G E, ie- rn fi n (x) = i/~ x G M*, n = l,...,N. (2.26) 

Our third requirement reads 

V!,..., v N G (0,oo). (2.27) 



There are explicit examples available where the second assumption ( 2,26] ) and the assump 



tion of Lemma ^lj are satisfied, but the third assumption ( 2.27 ) and symmetry of A are 
violated. These examples can be gleaned from Ref. ||, but we do not spell out the details 
here. 

Our fourth and final restriction can be most easily phrased in terms of the r-function 

t(x) = |1jv + CD(xY 1 \ . (2.28) 
In view of our second restriction ( |2.26| ), t(x) is an entire function, cf. (|2.6| )-( ^?7[ ). Likewise, 



from (2^) we see that R{x) can be written 

R(x) = E(x)/t(x), (2.29) 

where E(x) is an entire function. Our fourth requirement is now that t(x) have no zeros 
for Imx 6 [—1,0]. This entails in particular that R(x) has no real poles (the assumption 
of Lemma |2.1| ) . 

In Appendix A we prove that when N + or A_ equals N (so that all r n lie either on the 
positive or on the negative imaginary axis), the third restriction entails the fourth one, 



cf. Lemma 4.1 . For A + A_ > this is presumably still true for generic spectral data, but 
we were unable to prove this. Explicit examples we do not present here show that our 
fourth requirement is stronger than the third one. In any event, in Appendix A we also 



prove that for A + A_ > the third restriction (2.27) together with the requirement 

-ir x , -ir N+ ,-ir N++ i + tt, . . . , -ir N + n G (0,7r/2] (2.30) 
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are sufficient to obtain the fourth one, cf. Lemma A. 2 



We continue by showing that our fourth assumption suffices for essential self-adjointness 
of A on its definition domain V . In the proof we use one property of t(x) that cannot be 
found in Parts I and II, namely, 

t*(x) =t(x -i). (2.31) 
This formula is an easy consequence of ( 2,28| ) and the relations 

C = -C\ D*(x) = -D(x-i), (2.32) 
which follow from ( p.26| ). (Cf. also Appendix A and I Appendix D.) 

Lemma 2.2. Assume that the data satisfy ( \2.2fl( )~( 2.27 ) and that t(x) ( \2.2Sj) has no 
zeros for Imx E [—1,0]. Then the operator A defined by ([l.2£ ) is essentially self-adjoint 
on V \1.2\ ). 

Proof. As already detailed, it suffices to prove symmetry of A on V . (Recall the paragraph 
containing ( 2.25[) .) For this purpose we fix <fii,<p2 £ Co°(^) an d consider [T<f>\ : ATfa). 
By definition, this equals 

(FfaFMfa) = -!- lim f dx( f dqW(x,q)M<l) 
2tt R^-co J_ r 



dpW(x,p)2 cosh (p) 02 (p) 

oo 

oo r R 



= — lim / dq<})i(q) / dp<j) 2 (p) dxW(x,q)2cosh(p)W(x,p), (2.33) 
2itR^odJ_ od J_ od J_ r 

where we used Fubini's theorem in the last step. Rewriting (AFfa, J- fa) in the same way, 
we obtain 

{Tfa,ATfa)-{ATfa,Tfa) = — lim / &q<f> x {q) \ dp<p 2 {p)I R {p,q), (2.34) 
where 

f R 

Ir(p,0)= dxW{x,q)W(x,p)(2cosh(p) - 2cosh(< ? )), p,qeR. (2.35) 
J-R 

Next, we invoke the eigenvalue equation (|ll^) and the notation Ql.S| ) to rewrite Ir(p, q) 



as 



dx(W*(x, q)[W(x -i,p) + V a (x)W(x + i,p) + V b (x)W(x,p)] 

R 

- [W*(x + i,q) + V*(x)W*(x-i,q) + V b *(x)W*(x,q)]W(x,p)). (2.36) 
Recalling ([L~7|), we obtain 

f R 

Ir(p,q)= dx{J{x,p,q) - J(x + i,p,q)), (2.37) 

J-R 
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where 

J(x,p,q) = W*{x, q)W(x -*,p) - V a (x - %)W{x - i,q)W(x,p). (2.38) 
Now from ([Hi]) and ( gjg ) we have 

W(s,p) = f 1 - £ K{X) 1 - ) , (2.39) 



with E n (x) entire. Using ( [2.31 ), we deduce 



W*(x,q) = e- ix Ul-J2 f niX \ Q 1 r ), (2-40) 



From 11(2.34) we also have the identity 



t(x + i)r(x — i) 
r{xf 



V a (x) = '^'y^ . (2.41) 



When we substitute (|2T39|) (pT4l|) in J(x,p,q), we can write the result as 



N ™, > \ / N 



J(x,p, q) = e-™« 1-Y EUX \ 1 1 - V ^fe ~ f 

v y I r(x - «) ei -e r " J \ ^ t(x - i) eP - e~ r k 



1 / - 1 

. e -i(*-i)« r(x _2i)_^^( a 



x — z 



r(x — i) 2 \ v ' ^-^ " v ' e 9 — e Tn 

v ' \ n=l 



The point of doing so is that this representation shows that J(x,p,q) has no poles on and 
inside the rectangular contour T in the x-plane connecting —R, R, R + i, —R + i. (Indeed, 
by assumption t{x) is zero-free for Imx £ [—1,0].) 
As a consequence, the contour integral 

C R (p,q) = £dx,J(x,p,q) (2.43) 

vanishes by Cauchy's theorem. On the other hand, Cn{p,q) equals In(p,q) plus the 
integrals over the vertical sides of T. Thus we infer 

/ rR+i r-R \ 
iRfaq) =( + )dx 

\JR J-R+i/ 

x [-W*(x,q)W(x -i,p) + V a (x - i)W*(x - i,q)W(x,p)]. (2.44) 

In order to handle the right boundary term for R — » oo, we substitute 

W(x,p) = e ixp + p + (x,p), W*(x,q) = e- ixq + p* + (x,q), (2.45) 
V a (x) = l + p(x), (2.46) 
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in the first integral of ( |2.44| ). An inspection of (the proofs of) Lemmas 2A, 2J2 and 
Theorem 2.3 in Part I now reveals 

P+{x,p),p* + (x,q),p(x) 0, Rex^oo, (2.47) 

uniformly for p, q and Imx in M-compacts. An easy dominated convergence argument 
then shows that the contribution to ( |2.34| ) of terms containing at least one p vanishes. 
Thus we are left with 



R+i 
R 



_^—ixq+i(x—i)p _|_ ^ — i(x— i)q+ixp 



= -tejVp-*) sinh((p + g)/2) Sinh((p ~ g)/2) . (2.48) 

p-q 

Substituting this in t&M) , 

we can transform to sum and difference variables to infer that 
the contribution of this term vanishes by virtue of the Riemann-Lebesgue lemma. 
Next, we substitute 

W(x,p) = a(p)e ix 'P + p„(x,p), W*(x,q)=a*(q)e~ i ^ + p*_(x,q), (2.49) 

and (|2.46| ) in the second integral of (|2.44| ). Using the alternative representation 1(2.49) 
of a(p), we obtain 

N 



p-(x,p) = £ {{C-\) n - R n (x)) I (eP - e~ r ") 
n=l 
N 

= ]T (C- 1 D(x)R(x))J (eP - e"-) , (2.50) 

n=l 

where we used in the second step. As before, we now infer 

P-(x,p), p*_ (x, q), p{x) — * 0, Rex — > — oo, (2-51) 

uniformly for p, q and Imx in M-compacts. Thus it remains to consider the contribution 
of 

Aia*(q)a(p) exp[-iR(p - q)] sinh((p + q)/2) sinh((p - q)/2)/{p - q) (2.52) 

to ( |2.34| ). As before, this vanishes by the Riemann-Lebesgue lemma. Hence the rhs 
of ( p.34| ) vanishes, entailing symmetry of A on P. ■ 

With the assumptions of the lemma in effect, we can take the closure of the operator A 
on V to obtain a self-adjoint operator. The latter acts on a dense subspace of the Hilbert 
space 

Hx{F) = Ran(J c '). (2.53) 

It follows from the isometry of J-, which we prove in the next section, that the range of J- 
is actually closed. At this stage, however, we only know T is bounded and we have no 
information about Ran(^") and its orthogonal complement. 

Until further notice, we denote by A the self-adjoint operator on TC X that acts as the 
closure of A on Ji x {T\ and as an arbitrarily chosen self-adjoint operator on the orthogonal 
complement 7i x (J-)' L . The results of the next section are independent of the latter choice. 
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3 Time-dependent scattering theory 



Throughout this section, the assumptions of Lemma 2.2 are in effect. With the self-adjoint 
operator A on TL X defined at the end of the previous section, we show that the wave 
operators for the pair of dynamics exp(— itA) and exp(— UAq) exist and are intimately 
related to the eigenfunction transform T '■ As a corollary, this yields isometry of J- . The 
following lemma is the key to these results. 



Lemma 3.1. For all <ft E Cg°((0, oo)) we have 



lim || (T- J r o )exp(-itM)0|| = 0, 



t— >oo 

lim 

t— >— CO 



\(F - T Q a(-))e^{-itM)(j)\\ = 0, 



(3.1) 
(3.2) 



where M is given by ( LIS ) and a(-) is the operator of multiplication by a(p) \2.1h\ ). For 
all 4> G Cg°((— oo, 0)) we have 



lim - F a(-)) exp(-itM)<t>\\ = 0, 



t— iOO 



lim ||(^- F ) exp(-itM)<f>\\ = 0. 



(3.3) 
(3.4) 



Proof. To prove ( |3.1|) , we fix (j)(p) with supp(^) C [r, R], < r < R. Then we have 

a ixp— 2it coshp 



i roo N pR 

\\(F-F )exp(-ttM)4>\\ 2 = — dxY,Rn{x) dp 

Z7T J ~°° n =l Jr 



eP 



-4>(p) 



(3.5) 



Now when we change variables p — > y = coshp, we see that the integrals yield bounded 
functions b n (t, x) that converge to as t — > oo by virtue of the Riemann-Lebesgue lemma. 
By dominated convergence, this entails that the ^-integral over a bounded region converges 
to for t — > oo. 

To exploit this, we write 



dx 



dx+ dx + 



-i 



dx. 



(3.6) 



As we have just shown, the middle integral tends to for t — > oo. To handle the right inte- 
gral, we recall from Part I that R n (x) has exponential decay for x — > oo, cf. 1(2.41), (2.42). 
This decay supplies the domination we need, in combination with the pointwise conver- 
gence to of the x-integrand, to deduce it tends to for t — > oo, too. 

To handle the left integral, we use a stationary phase argument. Specifically, we write 



exp(ixp — 2itcoshp) = (ix — 2itsinhp) 1 d p ex.p(ixp — 2itcoshp), 



(3.7) 



and integrate by parts to get 

-1 N rR 



— [ dx ^^(x) f dpe lxp - 2ltcoshp d p 



1 



4>(p) 



[x — 2t sinhp) [eP — e 



(3.8) 
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Now we use the estimate 

\x - 2isinhp| 2 > x 2 + 4i 2 sinh 2 r, x G (-00, -1], pe[r,R], t > 0, 
and boundedness of R n {x) on R to obtain an upper bound of the form 

1 



C 



oo dX x 2 + ci 2 ' 



C,c> 0. 



(3.9) 



(3.10) 



The integrand is bounded above by the L ((— 00, — l])-function 1/x 2 and tends to as 
t — ► 00. Hence Q3.10 ) tends to as well, so that (|3.8[) does, too. Therefore, we have now 
proved (p.l|). 

In order to prove (|jO|), we observe that the estimate (|3.9|) is also valid for p G [— R, —r] 
and i < 0. Thus we can choose 4>{p) with supp(0) C [-R, —r] and proceed in the same 
way as for (|3.l| ). 

Next, we prove (U). To this end we recall ( gjgp and (H3>' which 

we rewrite as 



N 



W(x,p) - e ix Pa(p) = Y J Ux)/ {e p ~ e" r ") 



n=l 



with £ n (x) admitting the two representations 

e„(x) = (c- 1 c) n - Rn{x), 

e„(x) = (C- 1 J D(x)i?(x)) n . 
Now for (p(p) with supp(^) C [r, i?] we have 

||(J--^oa(-))exp(-itM)0|| 2 



(3.11) 



(3.12) 
(3.13) 



1 

2^ 



dx 



N 



n=l 



dp- 



ixp— 2it coshp 



-<t>(p) 



(3.14) 



From ( |3.12| ) we see that £ n (%) is bounded on R. Using once more the splitting (3.6), it 
follows that the middle integral tends to as t — > — 00. The right integral can be handled 
by the same stationary phase argument as before, noting that the estimate ( |3.9[ ) is also 
valid for x G [1, 00) and t < 0. 

It remains to show that the left integral tends to as t — > — 00. To this end we use the 
second representation (|3.13|) of £ n (x). Indeed, it entails that i n {x) has exponential decay 
as x — > — 00. (Recall the definitions fl2.6p , ( |2.7| ).) Thus we can once again combine the 
Riemann-Lebesgue lemma and the dominated convergence theorem to deduce convergence 
to for t — ► —00. 

Finally, to prove ( |3.3j ) we can proceed in the same way as for (3.2), noting that (3.9) 
also holds for x G [1, 00) and p G [-R, —r], M 

We are now in the position to obtain the principal result of this section. 

Theorem 3.2. The eigenfunction transform T ( 1-lSJ is isometric. The strong limits of 
the operator family exp(it^4) exp(— UAq)^ for t — > ±cxd exist and are given by 



U± = FA±{- 



(3.15) 
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where 

A+(p) = 
A-( P ) = 

The S-operator 

equals the unitary multiplication operator 

T{p) 

Proof. We choose ip(p) G C§°(R*), so that 



1, P>0, 

l/a(p), p<0, 

l/a(p), p>0, 

1, p<0. 



S n = U1U- 



l/o(p), p>0, 
a(p), p<0. 



1> = </>+ + 



0+ G C °°((0, oo)), 0_ G C °°((-oo, 0)). 



(3.16) 
(3.17) 

(3.18) 
(3.19) 

(3.20) 



By virtue of ([Sj) and ( |3.3|) we have 

\\Fip - exp(iti) exp(-iti o )-^o(0+ + a(-)<P-)\\ < Wi? ~ exp(z'iA) exp(-^i )-^o) 
+ ||(J r -exp(iti)exp(-iti )J r oa(-))</'-|| < {{(F - Fa) exp(-itAf)0+|| 
+ ||(jr-JF o a(-))exp(-itM)0_ 

From this we deduce 



0, 



i — > oo. 



(3.21) 



W^ipW = lim \\exp(itA)exp(—itAo)J : o( 

t— >oo 

= ll^+ + a(-)<M = II 



+ a 



(3.22) 



Since C^°(1R*) is dense in Ti p , it follows that T is an isometry. From (3.21) and its analog 
for t — > — oo we also obtain the second assertion of the theorem. The last assertion is then 
clear from ( |3.15j ) and isometry of T. M 

4 Bound states and spectral resolution 

We begin this section by focusing on the j4-eigenfunctions W(x,r k + 2ml), with x G R, 
k = 1, . . . , N and I G Z, assuming only (p])- (|2ll) . Using (glU) and (^5|)-(gl8|) we see 
they can be written 



W(x, rfc + 2ml) = exp(ir k x — 2irlx)d(rk, fj,p.; x)R k (x) 
Now from ([2.8D we readily obtain 



Z G 



lim d(r k ,fi k ;x)R k (x) = 1, fc = 1, ...,JV. 

Rex— >oo 



(4.1) 



(4.2) 



For square-integrability of ( [4.1| ) near oo we should therefore take I G N when Imr^, G (0, ir) 
and Z G N* when Imr fc G (— m 0). 
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Consider now square- integrability near — oo. Since Rk{x) tends to (C _1 C) fc G C* for 
x — > — oo (cf. ( |2.12D and I Lemma 2.1), we see from the definition (|2.7|) of d that we need 
—I € N. For Imrfc E (— 7T, 0), therefore, we cannot simultaneously have square-integrability 



of (4.1) near oo and near — oo. In contrast, for Imr^ £ (0,vr), the choice I = ensures 
square-integrability near ±oo. 

Thus far, we have not imposed restrictions on (r, fi) beyond our standing assumptions 



|2.1| )-(2.4). But to ensure square-integrability over R of W(x,r/%) for Imr/% € (0,7r), we 



should obviously require absence of poles for real x (the assumption made in Lemma 2.1). 
Doing so, we obtain j4-eigenfunctions whose restrictions to the real axis are in TC X . We can 
only prove pairwise orthogonality of these functions, however, when we make the same 
assumptions as in Lemma |2.2|. 



Lemma 4.1. With the assumption of Lemma \2. 1\ in force, the functions 

ip n (x) = W(x,r n ) = (j, n (x)exp(-ir n x)R n (x), n = l,...,N+, (4.3) 

satisfy 

VVi(^) ~ exp(ir n x), Rex — > oo, (4-4) 

ip n (x) ~ c n (C _1 C)nexp(-ir n x), Rex -> -oo, (4.5) 
uniformly for lmx in compacts, and their restrictions to R belong to 7i x . Now suppose 



that the requirements of Lemma 2.i are met. Then the functions ipi(x), . . . , ipN + {%) aw 
pairwise orthogonal in 7i x . 



Proof. We have already shown square-integrability and the asymptotics ( fl.4p -( |4~q ). To 
prove pairwise orthogonality, we invoke the eigenvalue equations 

(AW)(x,r n ) = 2cosh(r n )W(x,r n ), n=l,...,N+. (4.6) 

From ( |2.lD , ( |2.2D and ( |2.26[ ), we see that the eigenvalues are real and distinct. Letting 



n m, we now use (L6) to write 

2[coshr m - cosh r n ](V> n , ip m ) 

/oo 
dx(ip*(x)[ip m (x -i) + Va{x)ifj m (x + i) + V b (x)ip m (x)} 
-oo 

- [r n (x + v:(x)r n (x -%) + v b *(x)rjxM m (x)). (4.7) 

Using ( |1.7| ), we can write the rhs as 

dx(J nm (x) - J n m{x + i)), (4.8) 



with 

Jnm(x) = Ipni^mix - i) - V a (x - i)^* n { x ~ i)lpm(x). (4.9) 

We now recall ( |2.29| ). It entails we may write 

ip n (x) = fj, n exp(-ir n x)E n (x)/T(x), n = l,...,N + , (4.10) 
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with E n (x) entire. Using also ( 2.4l| ) and (|2.3l| ), we deduce that J n m{x) can be rewritten 

E^(x)E m (x - i) 



as 



Jnm(x) = V-n^m exp(-ir n x - ir m {x 



t(x 



exp(-ir n (x - i) - ir m x) 



E*(x - i)E m (x) 



(4.11) 



t(x — i) 2 

This representation shows that J nm (x) has no poles on and inside the contour T defined 



below ( |2.42| ) in the proof of Lemma |2.2| . Thus we can use the same reasoning as in 
that proof to infer that (4.8) vanishes. (The vanishing of the boundary terms is here a 
simple consequence of the asymptotics ( [4.4D , ( [4.5[) and Q1.2Q .) Hence pairwise orthogonality 
follows from vanishing of the lhs of ( [4.7D . I 

The lemma just proved together with our next lemma are the key to clarifying the 



character of the subspace Ran(jF) , when the assumptions of Lemma 2.2 are satisfied. 
From Theorem |3.2| we already know that in that case T is an isometry, so that we have 



(4.12) 



with P the projection on Ran(J 7 )^. In the next lemma we obtain a formula for (J 7 * fi, J-* f 2 ) 
with fx, f 2 € Cq°(R), from which this projection can be explicitly obtained. To prove the 
pertinent formula, however, we need only impose our second requirement ( p.26| ), together 
with absence of poles for real x, cf. Lemma [P|. 



Lemma 4.2. Assume that the spectral data satisfy ( 2.2t ), and assume R(x) has no real 
poles. Then we have for all fi, f% E Cq°(M) 



N+ 



(FfuFft) = (/i,/ 2 ) -^^(/l,^)(V>n,/2), 



(4.13) 



n=l 



where ^ n is defined by 

Proof. Our starting point is the formula 



j roc poo 

(r*f 1 ,F f f 2 ) = — lim / dyh{y) dxf 2 (x)l R (x,y), (4.14) 



where 



R 



R 



d P W*(x,p)W{y,p), 



x,y G 



(4.15) 



We are going to exploit that W(x,p) is the product of the plane wave exp(ixp) and 
a function of p that is meromorphic and 27rz-periodic, cf. ( p. 11 ). 

For this purpose we define the rectangular contour C connecting —R, R, R + 2iri and 
—R + 2iri in the p-plane, and put 



lc{x,y) = <j> dpW*(x,p)W(y,p), 



x,y G 



(4.16) 
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(Of course, the * refers here to the p-dependence, cf. fll.q) .) Due to our requirements on 
ri,. . . ,rjv, the only singularities of the integrand inside C consist of simple poles at the 
2N distinct points 

p = r n ,2-Ki - r n , n = l,...,N + , (4.17) 

p = -r n ,2iri + r n , n = N + + 1,...,N, (4.18) 

on the imaginary axis. Hence Cauchy's theorem yields 

l c (x,y) =27riTZ(x,y), (4.19) 

where 7Z(x, y) denotes the sum of the residues. 
On the other hand, we can also write 

Tc{x, y) = 1r(x, y) - e 2 ^ x - y h R (x, y) + B R (x, y), (4.20) 

with 

/ i-R+2-wi i--R \ 

B R (x,y) = ( + )dpW*(x,p)W(y, P ). (4.21) 

\JR J-R+2mJ 

Hence we have 

-i 



l R {x,y) = (l- e 2 ^-y>) [2mK{x,y) - B R (x,y)}. (4.22) 
We proceed to calculate the contribution of the residue sum to the inner product 

{F*h,F*h) = — lim / dyMJ) dxf 2 (x) 

1 - e^'A^ [2mK(x,y) - B R (x,y)], (4.23) 



where we combined ( |4.14 ) and ( |4.22| ). First, we show that the residues at the poles ( [1.18 ) 
cancel pairwise. Indeed, fixing I € {N + + 1, . . . , N}, the residue sum at p = —n, 2tti + n 
of the pertinent function 



\ fc=l / \ m=l / 

equals 

+ eMdn - 2*)(y -)l(-5M)(l-E ■ (1.25) 



Recalling the system ([Ts|), we see that this equals 



-exp[-iri(y - x)]d*{r h m\ x)R* (x)e ri Ri(y) 

-exp[{ir l -2ir)(y-x)]e- ri Rt{x)d(r l , t i l ;y)R l (y). (4.26) 
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From the definition (|2.7| ) of d we see that this is proportional to 



exp[— iri(y — x)\m exp[— 2i(r\ + iir)x]e ri 
+ exp[(ir z - 2n)(y - x)]e~ n m exp[-2i(r; + in)y] 
= exp[-irj(y + x)]e 2nx (jTie ri + nie~ n ) , 



(4.27) 



which vanishes due to ( 2.26| ). 

We are therefore left with the residues of ( 4.24j ) at the points ( |4.17 ). The residue sum 
at p = r n , 2ni — r n for n € {1, . . . , N+} equals (using (^8|) , (|2.26| ) and (fQjj)) 



exp[ir n (y - x)} 



R* n (x) 



-2ir n y 



R n {y) 



+ exp[(-ir„ - 2n){y - x)]72^e -ji,r " a R* n (x 



R n (y) 



exp[-ir n (y + x)]R* n (x)R n {y) ^ n e~ r " + e 2 ^ x ~ y ^e rn 
= iu- 1 (l - e 2nix - y) ^j e' irnX R* n (x)e- ir " y R n (y) 
iuJl-e 2 ^ x ~^)jJx)My)- 



(4.28) 



Substituting this in ( [4.23 ) and comparing to ( 4.13j ), we deduce that it remains to prove 



— lim 

2n R^oo 



dyh(y) / dxf 2 (x) 



Bn(x,y) 

,2-n{x~y) _ i 



(/l,/ 2 ). 



(4.29) 



To this end we first rewrite Br(x, y) Q4.21|) as 



Bn(x,y) 



R+i 



R—iir 

n(x-y) 



dp[W*{x,p + in)W(y,p + in) - W*(x, -p + in)W(y, -p + in)} 

R+iw 



dp 



R—in 



Mv-*)A{p, x, y) - e ip( - x - y ^A(-p, x, y) 



(4.30) 



where the auxiliary function A is given by 

A(p,x,y) 
We now claim that the identity 



\ n / \ m 



R m (y) 



A(p, x, x) = A(—p, x, x) 



(4.31) 



(4.32) 



holds true. 

To prove ( [4.32[) , we observe that the functions A(±p,x,x) are 2-7ri-periodic in p and 
bounded for |Rep| — > oo, and that they have simple poles in the period strip Imp € [0, 2n} 
at p = in ± r n , n = 1, . . . , N. Using (2^8) and ( |2.7| ) in the same way as above (cf. (|4.24j )- 
( 4.27] )), we readily verify that the functions have equal residues. By Liouville's theorem, 
it now follows that A(p, x, x) — A(—p, x, x) does not depend on p. 
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To show that this difference vanishes, we need only compare the Rep — > oo limits of 
A(±p,x,x). Obviously, we have 



lim A(p, x, x) = 1, 

Rep^oo 



lim A(—p,x,x)=X*(x)X(x), 

Rcp^oo 



where 



A' 



A(x) = 1 + eTmR ^ 



m=l 



Now from I(D.17) we see that (2.26) entails 

A*(x)A(x) = 1. 

Hence our claim (|3|) follows. 

Next, we introduce functions B and C by setting 

A(V,X,V) = 1 + e~ p B(p,x,y), 
A(-p, x, y) = A*(x)A(y) + e~ p C(p, x, y). 



Due to ( 4.32 ), these functions are related by 

B(p,x,x) = C(p,x,x). 

We now rewrite ( |4.30| ) as 

B R (x,y) = (B% ) (x,y)+B%\x,y) 

with 



R+i-w 



R+tTT 



dp 



R—in 



Jp(y~x) _ e iK*-»)A*(x)A(y) 



R-i 



dpe- p \e ip(y -' x) B(p,x,y) - e ip{x ~ v) C{p,x,y) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 
(4.37) 



(4.38) 
(4.39) 

(4.40) 
(4.41) 



The point is that we can get rid of the remainder term B^{x, y) by using ( 4.3Sj ). 
Specifically, ( 4.38 ) entails we may write 

i ^ f-R+iir f-y , . 

B%\ x ,y) = dpe~ p / dsd s (e ip ^- x) B(p,x,s)-e ip{x - s) C(p,x,s)) . (4.42) 

Now from the definitions flPfiD , ( ^37|) of B and C we deduce that for all s- values between x 
and y we have 



q s (fM»-*)B(p,x,s) -e^ ffi - s )C(p,z,s))| < \p\D(x,y), 
Rep > R, |Imp| < it, 



(4.43) 



where D(x,y) is a positive function that is bounded for x, y varying over M-compacts. 
Thus we obtain 



B%\x,y) < 2vre- K {R 2 +vr 2 ) 1/z |y - x\D(x,y). 



(4.44) 
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Hence the contribution of B^\x,y) to the lhs of ( 4.29| ) vanishes. 
We are now reduced to showing 



/OO l>00 »W 

dyhiv) I dxf 2 (x) . f 
-oo J-oo smh(vr 



4tt(/i,/ 2 ). 



r(x - y)) 

Calculating the integral ( 4.40| ), we can write the result as 

Bf (x, y) = 2 sinh(vr(x - y))[C c (R; x, y) + C S (R; x, y)], 



where 



C c (R;x,y) 



1 - A 



i(y -x) J 
C s (R;x,y) = (1 + \*(x)X(y)) 



cos(y — x)R, 
sin(y — x)R 



y-x 



(4.45) 

(4.46) 

(4.47) 
(4.48) 



By virtue of ( [4.35 ) and the Riemann-Lebesgue lemma, the contribution of C c (R;x,y) 
to ( [4.45 ) vanishes. Using the tempered distribution limit 



smcx 
lim = TT0[X) 

c^oo x 



(4.49) 



and (|4.35|) , the remaining term C s (R;x,y) yields ( |4,45[ ). 
Clearly, we can rewrite ( |4.13 ) as the operator identity 



N+ 

JTJT* = 1 - Vnlpn (g> 1p n , 
n=l 



(4.50) 



(4.51) 



where the rank-one operator t/j (g> Xi with ip, x G is defined by 

(V>®x)/ = (x,/)V>, 

For the remainder of this section, we assume that the requirements of Lemma |2.2| are 
met. Then T is an isometry (as proved in Theorem ET2), and the functions tp\, . . ., i(^n + ( [4.3| ) 
are pairwise orthogonal in 7i x (as proved in Lemma [Lll) . In view of ( [4.50| ), the projection P 
on Ran(^ r )- L can be written 

N+ 

n=l 

In particular, this entails the norm formula 



[il>n,ipn) = 



n = 1,...,JV+ 



(4.53) 



We are now in the position to turn the provisional definition of the self-adjoint Hilbert 
space operator A (see the end of Section 2) into a final one: We define A on Ran(^ r ) ± by 
(linear extension of) 



Aip n = 2cosh(r n )ip n , 



n = 1, 



(4.54) 
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Theorem 4.3. The operator A is essentially self-adjoint on the dense subspace 

C = V ffi Span(V>i, . . .,ipN+), 



(4.55) 



and its action on C coincides with that of the A AO A. The operator A has absolutely 
continuous spectrum [2, oo) with multiplicity two, and point spectrum {2 cosh(ri), . . . , 
2cosh(r7v + )} with multiplicity one. 



Proof. The first assertion follows by combining Lemmas 2A_ and 2J2 with (|4.3| ), ( [4.6[) 
and (|4.54 ). Denoting the domain of A by T>, the restriction of A to T> n Ran(J r ) is 
unitarily equivalent to multiplication by 2 cosh (p) on Ti p , cf. (1.22). Together with the 
distinctness of the numbers n, . . . ,rjy + G i(0,7r), this entails the second assertion. ■ 

We conclude this section with some further observations concerning three special cases. 
Taking first iV_ = 0, we recall that the assumptions 



< —in < 



< —irjy < IT, 



v n G (0, oo), 



n 



1,...,N, 



(4.56) 



are sufficient for all of our results to hold true, cf. Lemma |A.l . Now in this special case 
the AAO A satisfies 



where S+ is the AAO 



S 4 



A = S 2 + -2, 

exp(-id x /2) + V{x) exp(id x /2), 
V(x) = A(x)/A(z + */2), 



with A(x) given by ( 4.34 ); moreover, 

(S+W)(x,p) 



3 P/2 + e -p/2^j W(x,p). 



(4.57) 

(4.58) 
(4.59) 

(4.60) 



(These assertions follow from I Theorem 3.3.) Using the above results, we can associate a 
self-adjoint operator S + to S + by setting 

S+Ff = FM + f, (M + /)(p) = 2coBh( P /2)/(p), / €V(M+), (4.61) 
S+Vn = 2cosh(r n /2)Vn, n = l,...,N+. (4.62) 

It follows just as for A that the dense subspace C (4.55) is a core for 5+ on which the 
5 + -action coincides with that of S+. 

Secondly, we consider the special case = 2M, N + = N_ = M, together with spectral 
data 

< —ir% < ■ ■ ■ < —iru < k/2, fj G (0, oo), 

rN-j+i = Tj - in, VN-j+i = Vj, (4.63) 

where j = 1, . . . ,M. Again, this suffices for all of the above Hilbert space results to be 
valid, cf. Lemma A. 2. In terms of the particle variables defined in Appendix A, this choice 

-</.i/-, + i- "m-j-i-i- J = 1,...,M. (4.64) 



of spectral data amounts to 
< q+ < • • • < q+, i 
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Its distinguishing feature consists in the potential Vb{x) being identically zero. Moreover, 
after taking M — > N and performing a scaling x,p — > 2x,p/2, the class of AAOs A, 
together with their reflectionless eigenfunctions and associated self- adjoint operators A, 
amounts to the class of AAOs S+, together with their reflectionless eigenfunctions and 
associated self-adjoint operators S+. Once more, this follows from I Theorem 3.3. (The 
ordering we used there is different, but this is inconsequential. Indeed, all of the pertinent 
quantities are permutation invariant.) 



Thirdly, we consider the special case N+ = 0. From Lemma A.l we then infer that the 
assumptions 

-it < -ir N < ■ ■ ■ < -in < 0, i/„e(0, oo), n = l,...,N, (4.65) 
suffice for the validity of our Hilbert space results. In this case we have the AAO identity 

A = + 2, (4.66) 

where 

S- = exp(-id x /2) - V(x) exp(id x /2), (4.67) 
and V(x) is again given by (|4.59| ); moreover, 

= (e p/2 - e~ p / 2 ) W(x,p). (4.68) 

(These assertions are also a consequence of I Theorem 3.3.) Since T is unitary when N + 
vanishes, we can define a self-adjoint operator SL by 

S-=FM„F*, (M_/)(p) = 2sinh(p/2)/(p), / € 2?(M_). (4.69) 

As before, the subspace C = V is a core for S-, on which the SL-action coincides with 
that of S-. 

We would like to point out that the absence of bound states for the case = is 
a quite remarkable feature. Indeed, for reflectionless self-adjoint Schrodinger and Jaco- 
bi operators, absence of bound states implies that the potentials are trivial (constant), 
whereas here one obtains an infinite-dimensional family of nontrivial potential pairs V a ,Vb- 
When one takes the time dependence introduced in Part II into account, this family of 
reflectionless self-adjoint operators without bound states yields the left-moving soliton 
solutions to the analytic version of the Toda lattice studied in Part II. 

We can also use the = special case to illustrate the ambiguity issue discussed in 



the introduction, cf. in particular the paragraph below ( 1.18 ). Let us begin by noting that 



all of the wave functions W(x,p) studied in this series of papers satisfy 

(DW) (x,p) = 2 cosh{2irx)W(x,p) , (4.70) 

where the dual AAO D is given by 

D = exp{2md p ) + exp(-2-irid p ) . (4.71) 

Indeed, this is plain from W(x,p) being the product of the factor exp(ixp) and a factor 
that is a rational function of exp(p), cf. ( [2.11 ). For iV + = the operator T is unitary, so 



we can define a self-adjoint operator D on H p by setting 

D = T*M X T, (M x f)(x) = 2cosh(2irx)f(x), feV(M x ). (4.72) 
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The action of D on the core J-* (Cq°(R)) now coincides with the action of the A AO D. 
(This follows in the same way as the analogous assertion in Lemma |2.2| .) 

The upshot is that we have associated to the free AAO D an infinite-dimensional 
family of distinct self-adjoint reflectionless operators D without bound states. (Indeed, 
the function X(x) (4.34) plays the same role for D as the function a(p) ( |2.13| ) plays for A.) 



Interchanging x and p and performing a scaling by 27r, we see that we obtain a similar 
family associated with the free AAO Aq ( 1.11] ), as announced. 



Finally, we would like to use D with N + > to illustrate that Hilbert space operators 
associated to AAOs may look symmetric at first sight, even when they are not symmetric. 
Indeed, the symmetry property is far more elusive than may be apparent from the above 
results. (For instance, symmetry is probably generically violated when the requirements 
of Lemma 12 . 21 are not met, cf. also our results in Ref. 101.) 



For this purpose we observe that whenever the assumption of Lemma [^J is satisfied, 
we may define an operator D on the subspace T* (Cq°(R)) via 

Df*f = F*M x f, f G Cq°(R). (4.73) 

(The point is that we have 

Ker(^*) n C£°(R) = {0}, (4.74) 
by the argument proving ( |2.15| ); hence D is well defined.) Now with the stronger require- 



ments of Lemma 2.2 in effect, T is isometric, so that T* (Cq°(R)) is dense in 7i p . For 
N + > 0, however, the densely defined operator D is not symmetric. 

This assertion can be verified in two ways, both of which are illuminating. First, we 



can proceed as in the proof of Lemma |2.2| to try and show symmetry. Doing so, we are 



led to investigate the residue sum of the function ( 4.24| ) in the strip Imp G [0, 2tt]. As we 



have seen below (4.24), this residue sum vanishes only when N + = 0, so that D is not 
symmetric for N + > 0. 

The second way in which symmetry violation for N + > can be established hinges 
on the finite-dimensionality of Ran(J r )^ already detailed above. Specifically, Ran(J r )- L is 
JV-f -dimensional, and this suffices to rule out symmetry for N + > 0. 

To see this, assume D is symmetric on T* (Cq°(R)). Then it follows as before that D 
is essentially self-adjoint on J 7 * (Cq°(R)), and also that we have 

exp(itD)T*f = T* exp{itM x )f, f G C °°(M), t G C. (4.75) 

Choosing t real and taking closures, this entails 

exp(itD)T* = T* exp(itM x ), (4.76) 

and so 

exp(itD) = T* fscp(itM x )F, t G R, (4.77) 

by isometry of T. From this we deduce that the unitary one-parameter group exp(itM x ) 
leaves Ran(^") invariant. Hence it leaves Ran(^ r ) _L invariant, too. Since the generator M x 
has solely continuous spectrum and Ran(J 7 )- L is iV + -dimensional, we must have = 0, 
as advertized. 
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Appendix A. Absence of r(x)-zeros for Imx £ [—1,0] 

In this appendix we state and prove two lemmas that have a bearing on the most restrictive 
(fourth) requirement made in Section 2. We recall that this requirement consists in the 
restrictions ( |2.26| )-( p.27| ) on the spectral data (r,fi(x)) given by the first paragraph of 
Section 2, and in the additional restriction that t(x) ( |2,28D have no zeros for Imx G [—1, 0]. 
Our first lemma shows that when iV+ or N- vanishes, there is no need for the latter 
restriction. 

Lemma A.l. Assume that the spectral data satisfy 

< -in < ■ ■ ■ < -ir N < 7r, u n e(0,oo), n = l,...,N, (A.l) 

or 

-it < -ir N < ■ ■ ■ < -in < 0, zv n e(0,oo), n = l,...,N, (A.2) 

where v n is given by ( $.2(\ ). Then r(x) does not vanish in the strip Imx £ [—1,0]. 

For N + N- > 0, however, it seems likely that zeros of r(x) in the critical strip are 
not excluded by (|2~2^)-(^2p. (In particular, for N+ = N- = 1, we can show that t(x) 
may have zeros at the strip boundaries.) But an extra restriction on r ensures absence of 
critical zeros. 

Lemma A.2. Assume that the spectral data satisfy N + N- > and that in addition to 
we have 

IT TT 

< -in <■■■ < -ir N+ < 2 > < - ir N < ■ ■ < -ir N+ +i < --■ (A.3) 

Then r(x) has no zeros for Imx € [—1,0]. 

To prove these lemmas we invoke Section 5 in Part II. To this end, we require from 
now on (ET2t|) -( fr27D . Then we put 

a t = ~ ir i> i = X ' • • • ' N +> a r = ~ ir N++i + TT, i = 1, . . • ,iV_, (A.4) 
so that 

a s n E (0,vr), n = l,...,N 5 , 5 = +,-. (A.5) 
Next, we define real numbers 

qf = ln(cot(a+/2)), j = l,...,N+, 

q~ = -ln(cot(a^/2)) 5 i = l,...,JV_, (A.6) 

and positive 'potentials' 

V/(q) = J] |coth[(g+- 9 +)/2]| 

l<k<N+,k^j 

x \{ |tanh[(#- g f)/2]|, j = l,...,iV + , (A.7) 

1<Z<AT_ 

*T(?) = II |coth[(g--g-)/2]| 

l<m<Af_,m^Z 

x [] |tanh[(gf -g+)/2]|, « = 1,...,JV_. (A.8) 

1<3<N+ 
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Finally, we introduce real numbers 9j, j = 1, . . . , N + , 9 l ,1 = 1,..., iV_, by writing v n as 



cosh((/+ 



exp(0+), j = l,...,iV + , (A.9) 
2VT(V) 

^jv + +f = / , exp(flj-), l = l,...,J\L. (A.10) 
cosn(^ 



With this reparametrization of the spectral data in effect, the tau-function t(x) depends 
on points in the phase space 

« = ... G R 2iV | g+ + <■•• < S + 

<--'<<h>4 ±<h> 3 = h---,N + , Z = l,...,iV_} (A.ll) 

of the IIrei(V = tt/2) system studied in Ref. 0]. The crux is now that it can be rewritten 
as 

t(x) = \1 N + L(x + i/2)\, (A.12) 

where 

L(x) = £(q+,q~,9f -2afx, ...,0„_ -2a~ N _x), (A. 13) 

and C may be viewed as the Lax matrix of this integrable iV-particle system. Specifically, 
we have (cf. II Section 5) 

C(q, 9) = C{q + , q-)V(q + ,q-,9 + , 9~), (A.14) 
where T> is the diagonal matrix 

V ee diag(exp(0+)F+(<7), . . . ,exp(0+ + )V+ (?), 

exp(0r)Ff(g), . . . ,exp(^_)V^_(g)), (A.15) 

and C is the Cauchy matrix 

C jk = l/cosh[(qf -q+)/2], (A.16) 

C N+ +i,N + +m = 1/ cosh[(<?~ - q~)/2], (A. 17) 

C N++ i,k = -i/sinh[(gf - q+)/2], (A.18) 

Cj,N++m = i/smh[(qf - q~)/2], (A. 19) 

with j,k = 1, . . . , N+, I, m = 1, . . . , N-. We are now prepared to prove the lemmas. 
Proofs of Lemmas A.lLjA.2 . We fix £ G R and note that we may write 

r(£ - i/2 + mj) = |ljv + £(6^)1, (A.20) 
where U {rj) is the diagonal unitary matrix 

U{rj) ee diag(exp(2ia^7i), . . . , exp(2za^_??)), 77 G R. (A. 21) 
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Putting 

F( V ) = \U( V )+L(0\, (A.22) 

we should show that F(rj) does not vanish for any rj G [—1/2, 1/2]. For this purpose we 
first note that 

V(0 = diag(exp(0+ - 2a+£)y+, . . • ,exp(^_ - 2aJ r J)V„J ( A - 23 ) 

is positive. Similarity transforming L(£) with the positive diagonal matrix 2?(£) 1//2 , we 
obtain the symmetrized Lax matrix 

l s (0 = t>(0 1/2 cv(0 1/2 - ( A - 24 ) 

Since U (rf) is diagonal, we have 

F(rj) = \U( V )+L s (0\- (A.25) 

Consider now the inner product (f,L s (^)f) for a nonzero f eC N . When N+ or N_ 
vanishes, the matrix C is positive, so that (f,L s (£)f) > 0. With ( |A.5| ) in force, we also 
have the inequalities 

81m(f,U(v)f)>0, 6 (0,1/2], / + 0, 5 = +,-. (A.26) 

Therefore, Im (/, (U (rj) + L s (£)) f) does not vanish for \rj\ G (0, 1/2] and / 7^ 0. This entails 
that U(i]) + L s (^) is regular for \rj\ G (0, 1/2]. Since 1at + L s (^) is also regular, F(rj) ( |A.25[) 
does not vanish for r] G [—1/2, 1/2], and so Lemma Al follows. 
The restrictions ( A.3|) in Lemma |A.2| amount to 

a*G(0,7r/2], n = l,...,N s , 5 = +,-. (A.27) 

Hence we deduce 

Re(f,U(v)f) >0, |r?|<l/2, / G C". (A.28) 

Now for N + N^ > it is no longer true that C is positive. But in view of ( |A.16[) -( |A.1S| ) 
we may write 

C = C++ + iC s , (A.29) 

with C++ + C positive and C s real and symmetric. From this it readily follows that we 

have 

Re(/,L s (e)/)>0, fjLO. (A.30) 



Combined with ( A. 2^ ), this entails that U(rj) + L s (£) is regular for |r/| < 1/2. Hence 
Lemma IA.2I follows. ■ 



We would like to add that this proof involves only superficial features of the sym- 
metrized Lax matrix. Presumably, the far more detailed information obtained in Ref. [0 
can be used to relax the requirement ( |A.3| ) . 
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